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This paper explores the effects that magnetic fields have on the viscous boundary layers (VBLs) 
that can form in neutron stars at the crust-core interface, and it investigates the VBL damping of 
the gravitational-radiation driven r-mode instability. Approximate solutions to the magnetohydro- 
dynamic equations valid in the VBL are found for ordinary-fluid neutron stars. It is shown that 
magnetic flelds above 10^ Gauss signiflcantly change the structure of the VBL, and that magnetic 
fields decrease the VBL damping time. Furthermore, VBL damping completely suppresses the r- 
mode instability for B > 10^^ Gauss. Thus, magnetic fields will profoundly affect the VBL damping 
of the r-mode instability in hot young pulsars (that are cool enough to have formed a solid crust). 
One can speculate that magnetic fields can affect the VBL damping of this instability in LMXBs 
and other cold old pulsars (if they have sufficiently large internal fields). 



I. INTRODUCTION 

Gravitational radiation emitted by the r-modes always 
tends to drive these modes unstable, while internal dis- 
sipation most likely completely suppresses this instablity 
in all stars except neutron stars (see Lindblom ||] 

for a review). Understanding this instability in neutron 
stars is currently of great interest because unstable r- 
modes could emit gravitational radiation with an inten- 
sity that is detectable by gravitational-wave observato- 
ries, such as the proposed enhanced version of LIGO [|. 
A necessary step in determining the detectability of grav- 
ity waves from the r-modes is to find the critical angu- 
lar velocity that triggers the onset of the instability as 
a function of temperature. This splits the problem into 
two temperature regimes: newborn hot young neutron 
stars [^|-^, and cold old neutron stars spun up by accre- 
tion, such as in low-mass x-ray binary systems (LMXBs) 
[^-^. The reason for this split is that below approxi- 
mately 10^ K, the interior of a neutron star is expected to 
contain regions of superfluid neutrons mixed with lower 
concentrations of superconducting protons and normal 
electrons [^0|-Q. For fluid models, bulk viscosity is the 
dominant form of dissipation for T > 10^ K ; while for 
T < 10^ K shear viscosity dominates, and mutual friction 
due to electron scattering off quantized neutron vortices 
probably has little effect on the r-mode s p^ . 

However, Bildsten and Ushomirsky |I6| showed that 
when a solid crust is present, the shear dissipation in the 
viscous boundary layer (VBL) that forms at the crust- 
core interface is by far the most important suppression 
mechanism of the r-mode instability. Their work was ex- 
tended by Andersson et al. , Rieutord |l8| , Levin and 
Ushomirsky |l9| , and Lindblom, Owen, and Ushomirsky 
pof . The state-of-the-art calculation now includes the 
effects of the Corelois force, angular structure, equation 
of state, and compressiblity on the VBL. However, mag- 



netic effects on the neutron star VBL have not yet been 
published. 

The purpose of this paper is to explore the effects 
of magnetic fields on the VBL in neutron stars and 
to investigate the VBL damping of the gravitational- 
radiation driven r-mode instability. Approximate so- 
lutions to the magnetohydrodynamic (MHD) equations 
valid in the VBL are found for ordinary-fluid neutron 
stars. Neutron stars are expected to form a solid crust 
for p < 1.5 X 10^'^g/cm^ and for tempertures below an 
approximate melting temperature of T = 10^" K |^ , p2| . 
Once the crust-core region cools below 10^ K superfluid 
effects could become important. However, uncertainties 
exist in the quantum nature of the fluids, and it is pos- 
sible either the protons, or the neutrons, or both fluids 
could be in the normal state at the location of the VBL 
(see Epstein for a review of some of these uncer- 
tainties). Thus, this paper will focus on the case of an 
ordinary-fluid core. The crust is treated as a perfectly 
rigid, infinitely conducting solid. 

In this paper, it is shown that magnetic fields above 
10^ Gauss significantly change the structure of the VBL, 
and that magnetic fields decrease the VBL damping time. 
Furthermore, VBL damping completely suppresses the 
r-mode instability for B > 10^^ Gauss. Magnetic field 
strengths of hot young pulsars are typically similar to 
this lower bound. Cold old pulsars, however, tend to 
have fields less than 10^"'^ Gauss. Thus, assuming inte- 
rior fields are of a similar magnitude to external fields, 
magnetic fields will profoundly affect the VBL damping 
of the r-mode instability in hot young pulsars (that are 
cool enough to have formed a solid crust). One can spec- 
ulate that magnetic fields can affect the VBL damping of 
this instability in LMXBs and other cold old pulsars (if 
they have sufficiently large internal fields). This possibil- 
ity, and how superfluidity might effect it, are discussed 
further in the conclusion. Also, the results of long term 
magnetic evolution on the interior field need to be con- 
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sidered when applying the results in this paper p3| , p4[ . 
Finally, the complications associated with nonlinear mag- 
netic effects, such as the azimuthal winding of field lines 
psi , should also be taken into account. Such effects are 
not likely to change the general conclusions presented 
here, though nonlinear effects are beyond the scope of 
this paper. 

The next section reviews the MHD equations for neu- 
tron stars. Sec. Ill finds an approximate solution to 
the MHD equations valid in the VBL. Sec. IV discusses 



limits on the magnetic field and the VBL length scales. 
Sec. ^ presents the results of numerical calculations for 
the VBL damping times and the critical angular veloc- 
ity for two models of the magnetic field. Checks on the 
self-consistency of the results are made at the ends of 
Sees. IV and M. Caveats about the approximations used 
to obtain the results are discussed, and comments on how 
this work could be extended, especially to include super- 
fluid effects, are made in Sec. VI. 



II. MHD EQUATIONS 

The MHD equations for a mixture of superfluid neu- 
trons, and ordinary-fluid protons and electrons were de- 
rived by Easson and Pethick More recently, the 
MHD equations for a mixture of superfluid neutrons, type 
II superconducting protons, and ordinary-fluid electrons 

based on the earlier work 
Mendell MM, and Lind- 



were derived by Mendell |27 
of Mendell and Lindblom [ 28 
blom and Mendell To study the ordinary-fluid case, 
the ordinary-fluid limit of the equations in Mendell ||2^] 
is taken |^^. In this limit, scattering between particles 
causes them to approximately flow as a single fluid. How- 
ever, the presence of magnetic fields does allow small but 
finite electrical currents to exist. Thus, the relative ve- 
locity between the protons and electrons is not exactly 
zero, but nearly so. These currents are accounted for by 
the MHD approximation. 

The ordinary-fluid limit of the equations is accom- 
plished by defining the average velocity, w, in terms of 
the mass densities and velocities of the neutron, Vn, and 
the average velocity of the charged particles (the protons 
and electrons), u, by 



pv 



+ PpU, 



(2.1) 



where p = pn + Pp, and pn and pp are the mass densities 
of the neutrons and protons respectively. For simplicity, 
terms of order p^/ Pp are ignored. The dynamic equation 
for V is found by taking the equations in Mendell |27| for 
v„ and u and performing this averaging. Furthermore, 
ignoring small corrections due to electrical currents, it is a 
good approximation to replace the electron velocity with 
the average velocity in the resulting equations. Finally, 
quantized vortices do not occur in ordinary-fluids, and 
the forces associated with them are not included. Thus, 
for a star rotating uniformly in equilibrium with angular 



velocity f2, the equation for v for small perturbations can 
be written in the corotating frame as 



P \ c p 



(2.2) 



The only other independent vector degree of freedom in 
the MHD limit is the magnetic field, B, given in the 
corotating frame by 



dt6B = V X ((Jwx B), 



V • (55 = 0. 



(2.3) 



(2.4) 



Equations (|2.2|)-( |2.4D are the standard textbook MHD 
equations, except the major contribution to the mass 
density comes from the neutrons. Only the dominant 
dissipative effect due to the shear of the fluid, (?, is in- 
cluded. When magnetic diffusion is important, magnetic 
effects lead to Hartmann type boundary layers . How- 
ever, note that Easson calculates that the ratio of the 
magnetic diffusion time scale to the viscous diffusion time 
scale in a typical neutron star is roughly 10^(10^" K/T)^. 
Thu s, th e magnetic diffusion term has been ignored in 
Eq. ( |2.3| ). This means the magnetic effects on the VBL 
found here are quite different from the Hartmann case. 
Finally, note that SU = Sp/p + S^, where p is the pres- 
sure and <I> is the Newtonian gravitational potential (see 
pl| and references therein). 

The MHD limit is valid for studies of oscillations with 
phase velocities much less than the speed of light, and 
frequencies much less than the plasma and cyclotron fre- 
quencies. Under these circumstances the above equa- 
tions, along with the mass conservation laws and equa- 
tions of state, completely determine the dynamics of the 
system. All the other vector fields of interest are de- 
termined in terms of 6v and SB. Specifically, for phase 
velocities much less than the speed of light, the displace- 
ment current can be ignored in Ampere's law, and the 
current density is given by 



SJ 



in 



-V X SB. 



(2.5) 



For an infinitely conducting crust, there will also be a 
surface current density at the crust-core interface, given 

by 
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c 

An 



6B X f 



(2.6) 



where Rc is the radius of the core. For finite constant 
conductivity, magnetic fields can diffuse into the crust 
via 



dtSB = — V^^B, 
47rl 



(2.7) 
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where T is the conductivity. However, while Eq. (2.7) 
is not too difficult solve, in this paper it will only be 
necessary to use this equation to place bounds on the 
electron drift velocity and the rate of ohniic dissipation 
in the crust. Returning to the discusion of the quantities 
in the core, the equilibrium charge density is negligible 
(as is the perturbed charge density for frequencies much 
less than the plasma frequency), and thus the number 
densities of the electrons and protons are equal: rig = 
np = pp/mp. Thus, the relative velocity between the 
protons and electrons is given by 



5vp — 5ve — — -SJ. 

epp 



(2.8) 



It will be shown near the end of Sec. IV that the frac- 
tional difference between the proton and electron velocity 
for the r-mode frequency is of the order 10~^° for typ- 
ical neutron numbers. This fact, along with the strong 
coupling between the protons and neutrons, justifies the 
approximation 



5v„ 



5v. 



(2.9) 



Finally, when the above approximations hold and the 
conductivity is high, electrons (being the least massive 
charge carrier) respond to make the Lorentz force on 
them negligible. Thus, again ignoring small corrections 
due to the current density, the electric field is given by 



5E = J-^y.B-^-y.5B. 
c c 



(2.10) 



III. 



APPROXIMATE VISCOUS BOUNDARY 
LAYER SOLUTIONS 



Approximate solutions to the MHD equations are 
found in this section. The equilibrium magnetic field 
is assumed to be arbitrary (for now) except that it is 
static in the corotating frame, and it is restricted such 
that no equilibrium electrical currents exist. This implies 
the equilibrium structure of the star is unchanged by the 
presence of the magnetic field. To facilitate the manipu- 
lation of tensor quantites, a rotating spherical coordinate 
basis will be used, and indices will be raised and lowered 
using the metric tensor. In this basis, the equilibrium 
velocity is — H.cj)"'. Note that, following the notation 
of previous studies, Latin indices are space indices, ex- 
cept where it is understood that n, p, e, and c refer to 
neutrons, protons, electrons, and the crust respectively. 

Let Sv'^ and SU describe the nondissipative, lowest or- 
der (when expanded in powers of O), r-mode solution 
ignoring magnetic fields. These solutions are valid in the 
bulk of the core where viscous and magnetic forces are 
small compared to the Coreolis force. However, because 
of viscosity, the no-slip boundary condition must be ap- 
plied to the velocity at the crust-core interface. This 



causes a VBL to form. In the VBL, the magnitudes of 
the viscous, magnetic, and Coreolis forces become com- 
parible. Thus, all these forces have effects on the struc- 
ture of the VBL. Let 6v'^ and SU describe the corrections 
that must be added to the standard r-mode solutions to 
enforce the no-slip boundary condition 



0. 



(3.1) 



Finally, note that for the approximation of an infinitely 
conducting crust, the only boundary condition that must 
be applied to the magnetic field is that it gives the surface 
current density given by Eq. (2.6). The surface current 
density is small, and does not significant ly c hange the 
no-slip boundary condition given by Eq. (3.1) (see near 
the end of Sec. |[V|) . In this case, no magnetic diffusion 
takes place, and the details of what happens in the crust 
can be ignored. 

It is possible to find approximate equations for the cor- 
rective quantities by taking the following steps. First, let 
all perturbed quantities in the corotating frame have time 
dependence exp(jKi7t). Second, note that since Sv"^ van- 
ishes everywhere, the boundary condition that the radial 
components of the velocity fields vanish at the crust-core 
boundary is satisfied if the radial component, ^C, also 
vanishes. Finally, following the work of Lindblom, Owen, 
and Ushomirsky |^^, and references therein, all angular 
derivatives of the corrective quantities are ignored. This 
is valid, since the no-slip boundary condition forces the 
corrective velocity to rapidly change in the radial direc- 
tion. Thus, terms involving radial derivatives dominate 
all others. (For certain angles, certain caveats must be 
added to this statement. These angles are discussed in 



Sees. IV and VL) Typically, the ratio of radial deriva- 
tives to angular derivatives is the ratio of the radius of 
the core to the boundary-layer thickness. Taking these 
steps, the resulting equations (expressed in a corotating 
spherical coordinate basis) are 



SB" 



in 



B'-drSv", 



(3.2) 



iKflSv" - 2f]cos0(sin05{i"^) = TdlSv\ 

where 



Vl (B, 



iflKB 



IT] 
P 



(3.3) 
(3.4) 
(3.5) 

(3.6) 



In Eqs. ( pT^ ) and (3.6) is the square of the Alfven 
wave speed, defined as 
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47rp' 



(3.7) 



and throughout this paper note that B — \B\. Note 
that the magnetic effects on the VBL vanish if the radial 
com pone nt of the equihbrium field, Br, vanishes. Also, 
Eq. (3^) for the quantities 5U is included for complete- 
ness, but plays no further role in this paper. Finally, 
Mendell has shown that Alfven waves are replaced 
by cyclotron- vortex waves in a superfluid neutron, type II 
superconducting proton, ordinary-fluid electron plasma. 
Thus, it is possible to speculate on what happens when 
quantum fluids are present. This will be done in Sec. [v| , 
where expected complications with the superconducting- 
superfluid case will also be discussed. 

To find solutions to the equations, first note that the 
(/)-component of the equations becomes identical to the 
(^-component when 



(3.8) 



The corrective boundary layer solution is then found 
by allowing all perturbative quantities to vary as 
exp[ifc(i?c — r)]. It is then easy to show that solutions 
exist for 



\ nil \ B J ^ p 



where 



K± = Vk ± 2cos6'. 



(3.9) 



(3.10) 



Choosing solutions where Im(fc±) > 0, so that the correc- 
tive quantities decay exponentially as r — > 0, the general 
solution for the ^-component of the velocity is 



(3.11) 



The 0-component of velocities follows from Eq. (3 



If we apply the no-slip boundary condition given in 
Eq. (3.1) the constants C± are given by 



(3.12) 



The components of the standard r-mode velocity in the 
corotating frame (with their time dependence cancelled) 
are |o| 



5v' 



-iAr'"-isin"-i0e™*, 



(3.13) 



(3.14) 



As shown in previous papers, the viscous damping rate 
is given by integrating the shear over the core: 



1 



1 

2E 



(3.15) 



where E is the energy of the mode as defined in e.g., 
Lindblom, Mendell, and Owen [|^, but limited to the 
core. The largest contribution to the integral comes from 
the radial derivatives of the corrective boundary layer 
velocity, so that 

Sal.Sa'^" = ^Rli\drSv'\^ + \dr{smeSi}^)\'), (3.16) 

ignoring terms smaller than these by a factor of the 
boundary layer thickness over the core radius. Thus, it 
can be shown that the solutions presented in this section 
give a VBL damping time of 



27r2'"+3(m + 1)! f"" ( r 



■ql m(2m + 1)!! Jq 



Rc 



2m+2 



dr, (3.17) 



where 
1 = 



/""[|fc+pd+(l-COS0)2 

+ |A:_|2d_(l +cos6l)2 
Using the definition 

2tt i 
A± d± 



sm 



2m-l 



(3.18) 



(3.19) 



the boundary layer thicknesses in Eq. ( 3.18| ) and bound- 
ary layer oscillation wavelengths are defined by 



27r 



Re(A:±)' 



d+ = 



Ini(fc-i 



(3.20) 



(3.21) 



The results given here then reduce to those given by Lind- 
blom, Owen, and Ushomirsky 12^ in the limit B 0. 



IV. BOUNDS ON THE MAGNETIC FIELD AND 
BOUNDARY LAYER LENGTH SCALES 

In this section, numerical estimates are made for the 
characteristic length scales that determine the properties 
of the VBL. A lower bound is placed on B, above which 
magnetic field effects are dominant. Towards these ends, 
let Br = B for this section. 

The magnetic effe cts on the VBL are determined by 
examining Eq. (3.9). In this equation, if the magnetic 
term is larger than the viscous term, then magnetic ef- 
fects on the VBL length scales will be important. Thus, 
the condition for this to happen is 



Kf2 



> 



V 



(4.1) 



4 



Defining pi.^eii = p/(l-5 x lO^^'g • cm Tw = 
T/{10^'^ K) (and so on for other numeric subscripts in cgs 
units thoughout the rest of this paper) the ordinary-fluid 
viscosity is given by 



77 = 2.73 X 10' 



9/4 rji 



10 



(4.2) 



Using Eqs. (3.7) in Eq. (4.1) yields the following lower 
bound on the magnetic field, such that magnetic effects 
dominate the VBL properties: 



T-l 

-10 ■ 



(4.3) 



Note that hot young ordinary-fluid neutron stars have 
external fields much larger than the lower bound given 
above. Thus, magnetic fields will profoundly affect the 
VBLs in these stars if the internal fields are of a similar 
magnitude. Furthermore, note that the range of values 
of B in LMXBs probably lie on either side of the lower 
bound given above. Thus, it is uncertain to what extent 
magnetic effects on the VBL are important in these stars 
or other cold old pulsars with weak magnetic fields. Fur- 
ther complications with interpreting the results presented 
here in cold old pulsars occur due to the likelyhood that 
the protons and neutrons condense into superconduct- 
ing and superfiuid states. It has already been pointed 
out that Alfven waves are replaced by cyclotron-vortex 
waves in that case 



|27 



More is said about how this 
might change the results in this paper for cold old pul- 
sars in Sec. Vl. 



The characteristic length scales of the VBL are deter- 
mined by examining Eq. (|3.9|) an d the definitions of the 
length scales given in Eqs. ( p.20| ) and ( |3.21 ). Dropping 
the "±" symbol for the purposes of this discussion, since 
the VBL damping rate is controlled by terms of the form 
ri\k\^d — ?7(47r/A^ -\-l/(P)d, it would seem that the small- 
est length scales would determine the VBL damping rate. 
However, the ratio d/X turns out to be more important 
when magnetic effects are important. This point will be 
elaborated on in the next section. Also, only real frequen- 
cies, kU,, will be considered. This is valid either when the 
system is driven at a real frequency, or when the imagi- 
nary part of the frequency is small and can be ignored, 
except near certain special angles. The latter situation 
is the case for the r-modes and most other neutron star 
oscillations of interest. In this case, the coefficient K± is 
either purely real or purely imaginary depending on the 
angle 9. The more complicated case that occurs near the 
special angles already mentioned is discussed at the end 
of this section and in Sec. Also, to further simplify 
the discussion, it will always be assumed that k > 0. 

Taylor expanding Eq. (^|^) for the case when B is less 
than the lower bound given in Eq. ( |4.3| ) , the viscous term 
dominates, and the boundary layer thicknesses and wave- 
lengths are given to lowest order by 



27r \K±\\ np 



(4.4) 



For typical neutron star numbers these length scales are 



2^ 



0.024 cm 5/, 



Pl!5el4'^10 ^20d07r 



1/2 



(4.5) 



These equations agree in form with the standard result, 
but differ from previous results by numerical factors of 
order unity that are incorporated into \K±\. 

Taylor expanding Eq. (B.£) for the ca se w hen B is 
greater than the lower bound given in Eq. (L3), the mag- 
netic terms dominate, and the wave number is given by 



k+ = K. 




i r] K^l 



(4.6) 



The problem now divides in two cases corresponding to 
when K± is real and imaginary respectively. 

For the case of real K±^ the boundary layer length 
scales are given to lowest order by 



As 
27r 



1 



\k±\k^/^ n ' 



(4.7) 



(4.8) 



The subscripts indicate whether these quanties depend 
on purely magnetic, or a ratio of magnetic to viscous 
quantities. Substituting in values for the parameters 
gives 



\b _ 0.017cm ^ 1 _i/2 

27r ~ |;^^|^l/2-'^4.6e9"20007rPl.5el4' 



di}/„ = 



0.033 cm 



B/ri ~ |j^^|^3/2-^4.6e9^'20007rPl. 5614^10 



-11/4^2 
J 1 1 



(4.9) 



(4.10) 



Obviously As and dg^^ d o no t appear much different 
from what is given in Eq. (4^) for the limiting field of 
B = 4.6 X 10^ Gauss. However, for a typical 10^^ Gauss 
field 



Ab 

27r 



3.7 cm 

|X±|k1/2 



£'12"20007rPl.5el4J 



(4.11) 



_3.4xl05cm 



-Sl2^20007i-Pl.5el4^U)- (4-12) 



Note that in Eq. ( [1.12 ) the boundary layer thickness is 
approaching 30% the size of the star. Thus, the approx- 
imation that d/Rc is small becomes less valid as the size 
of B increases. For the case of imaginary K±, the roles 
of A and d become interchanged. Thus, for t his case, the 
results are the same as for Eqs. (4.7) - (4.12) with 



A 2TTd. 



(4.13) 
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Two checks on the self-consistency of the solution can 
now be made. First, it can be shown that for the length 
scales presented here that the fractional difference be- 
tween the proton and electron velocities is neglibib le in 
the VBL, as has been assmned. Substituting Eq. (3.2) 
into Eq. (2^) and the result into Eq. (2^) gives 



6vl - SiC = i(^ 
^ \ iripC 



( bv^ 



(4.14) 



which for typical neutron star numbers gives 
[(1.3 X 10-10 cm2)fc2]Bi2 



Sv" 



El) ( _ ^^^QSiH- 

B I \si\yO 



(4.15) 



Thus, the fractional difference between these velocities 
is small as long as (10-^° cm^)A;^ <C 1, which is true 
for the length scales presented in this section. Second, 
at the boundary the electrons can sl ip t o produce the 
surface current density given by Eq. (2.6). The surface 



current will actually be spread out over a thickness equal 
to the skin-depth, that t he m agnetic field penetrates 
the crust. Examining Eq. (2.7), it can be seen that for 
the r-modes that C, is given by 



^ V 2TTnnV' 
The conductivity, F, is given by [ p6| , 

F = 2x lO^Ss-Vpisrfo', 



and thus. 



_ 3.4 x 10 ^cm ^_i/2 -3/4 

"t. - -f^ i'20007rPpl3 -'lO- 



(4.16) 



(4.17) 



(4.18) 



Within the skin-depth region of the crust 5v'^ — 
-{mp/epp)SJ'' and SJ" ^ SI^/C !§. Thus, the 
electron contribution to the average velocity at r = 
Rc, which is {pe/p)Sv^, is negligible as long as 
(10-10 cm2)(pe/p)(fc/C) < 1, which also holds true for 
the length scales presented here. Ohmic dissipation 
within the skin-depth region is considered at the end of 
the next section, once the dominant VBL length scales 
are determined. 

Finally, consider the special case that occurs near cer- 
tain special angles, defined by 



K ± 2cosX± = 0. 



(4.19) 



For regions of the crust-core boundary where 6 is near 
X±, the small imaginary part of k becomes impor- 
tant. To understand this, make the replacement k — > 



K + i/(f2r„), and note that flTy ^ 1 for the situations of 
interest. In this case, for 9 = X± the angular derivative 
dk±/d9 becomes proportional to y/ilTy/i, which is large 
(and infinite at X± when the imaginary part of k is ig- 
nored). However, these regions make little contribution 
to the VBL damping rate, since oc l/y/Hrv is small 
near X± . Thus corrections to the solutions near 9 ~ X± 
are not made in this paper, nor have they been made in 
previous studies. 



V. BOUNDARY LAYER DAMPING TIMES AND 
THE CRITICAL ANGULAR VELOCITY 

As shown in previous studies, gravitational radiation 
emitted by the r-modes always tends to drive these modes 
unstable. Of primary interest to this study is how mag- 
netic fields change the onset of the instability of the r- 
modes via their effect on the VBL damping rate. Since 
the VBL damping rate dominates all other forms of 
dissipation, the onset occurs when this rate equals the 
gravitational-radiation growth rate. Thus, we can cal- 
culate the critical angular velocity for the onset of the 
gravitational-radiation instability by solving 



(5.1) 



where Ty is given by Eq. (3.17). It has been shown that 
TGR for the r-modes is given by 



TQR = TGR 



n 



2m+2 



(5.2) 



where flo = ^/irGp, and p is the average density of the 
star. 

Attention is now restricted to the case used in pre- 
vious studies: the m = 2 r-modc for a 1.4Mq n = 1 
polytrope, as described in, e.g., Lindblom, Owen, and 
Mendell and Lindblom and Mendell [|l5|. Adoption 
of this model allows easy comparison with previous stud- 
ies. (Note also that the m = 2 case is the most suscep- 
tible to the r-mode instability.) For this case, kq = 2/3, 
the stellar radius is 12.53km, and flo = 8413 s-^. Note 
that the maximum angular velocity, where mass shedding 
occurs at the equator, is roughly 2r2o/3 (for any equa- 
tion of state). The density at the crust-core boundary is 
given approximately by 1.5x IQi^ g/cm^ (see pl| , p^ ). Us- 
ing this density, the characteristic gravitational-radiation 
growth time is tgr = 4.25 s [^, and the cor e ra dius is 
Rr = 11.01 km. The viscosity is given in Eq. (4.2). 



Two models of the magnetic field will be considered. 
For both models the overall magnitude of the equilib- 
rium magnetic field, B, is taken to be a constant in the 
boundary layer. The models are defined by 

Model I Bj. = B = constant. 

Model II Br = ScosO, (5.3) 
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TABLE I. The VLB damping time in seconds for fl = Qo- 
Model II numbers are for A — -k 12. 



Model I 
T = 10** K 



Model I 
r= 10" K 



Model II 
T = 10** K 



Model II 
T = 10" K 



0.0 

10^ 
10"' 
10^1 
10^2 



51.89 
51.89 
51.88 
50.57 
16.03 



5189 
5056 
1603 
161.0 
16.10 



51.89 
51.89 
51.89 
51.41 
27.06 



5189 
5141 
2706 
303.4 
30.37 



where Q is the angle between the magnetic axis and the 
vector to the point (6', </>). This is given in terms of the 
location of this point by 



cos8 — sinAsin^cosf/) + cosAcosfl, 



(5.4) 



where the magnetic axis is chosen to lie in the xz plane of 
the corotating frame, and A is the angle between this axis 
and the rotation axis. Note that the 9 and (j) components 
of the magnetic field cannot be zero for Model II if B is 
constant. However, these components do not enter into 
the effects on the VBL. 

Now note that the VBL damping time, r^,, is a func- 
tion of the temperature, angular velocity, and magnetic 
field, but in general the depend ence on these quantities 
cannot be factored out of Eq. ( |3.17 ) as they have been 
in previous studies (but see towards the end of this sec- 
tion). Still, one can numerically evaluate for various 
temperatures and models of the magnetic field. The re- 
sults are presented in Table | for O = fio. Note that the 
values of t„ agree with the results of Lindblom, Owen, 
and Ushomirsky in the _B = case. 

The temperature dependence of the critical angular ve- 
locity is presented in Figs. In these figures, the hor- 
izontal dashed line corresponds to 2f2o/3, which is the 
approximate maximum angular velocity for which mass- 
shedding occurs. Thus, the regions above this line are 
unphysical and are kept only to illustrate the dependence 
of the curves for the range of magnetic field magnitudes 
typically found in neutron stars. Also, the curves prob- 
ably only apply for T > 10^ K, since only ordinary-fluid 
neutron stars are considered in this paper. They have 
been extended below this value to make comparison with 
other studies easier. This also makes allowances for un- 
certainties in the superfluid transition temperature at the 
crust-core boundary. 

Figure^ shows the temperature dependence of the crit- 
ical angular velocities for Model I. Figure 2 shows a simi- 
lar result for Model II when A = tt/2. Figure 3 compares 
the B = case with the Model I and Model II (A = 
and tt/2) B = 10^^ Gauss cases. (Results for other val- 
ues of A fall between the two Model II curves shown.) As 
would be expected, the angular dependence of the mag- 
netic field for Model II reduces the overall effect of the 
magnetic field on the critical angular velocity, as com- 
pared with Model I. Overall, it is seen that for magnetic 




FIG. 1. Temperature dependence of the critical angular ve- 
locities for Model I and i? = 0, 10® (indistinguishable from 
B = 0), 10", 10", and 10^^ Gauss, from bottom to top. 




FIG. 2. Temperature dependence of the critical angular ve- 
locities for Model II and -8 = 0, 10® (indistinguishable from 
B = 0), 10", 10", and lO" Gauss, from bottom to top. 
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FIG. 3. Temperature dependence of the critical angular ve- 
locities for Model I B = Gauss (bottom curve) , Model I 
B = 10" Gauss (top curve), and Model II B = 10" Gauss 
for A = and A = 7r/2 (middle two curves). 




FIG. 4. Integrand that gives X for </> = 0, T = 10* K, 
Q. = Q,o, and B = (bottom curve) and Model I 
B — 10^^ Gauss (top curve). 



fields B > 10"'^^ Gauss that the r-mode instability is com- 
pletely suppressed in ordinary-fluid neutron stars. 

Note that the critical angular velocity becomes tem- 
perature independent for high temperatures. This can 
be understood by considering the length scales from the 
last section. It is apparent that for large B that the VBL 
damping rate varies as 



1 



which implies for large B that scales as 



le.lOsB 



12 ' 



(5.5) 



(5.6) 



of 



independent 
Eqs. (O), (p|)-(pO|)] 



temperature or angular velocity [see 
The constant of proportional- 



ity was found by fitting the data in the table for Model 
I. Substituting this equation into Eq. (^) and using 
Eq. (|]^) gives the following results for Model I critical 
angular velocities 



(5.7) 



The Model II result is identical to this, but with a co- 
effient of 0.72. This equation gives the exact numerical 
results for the regions where the critical angular velocity 
curves flatten out. 

Thus, the temperature independence of the critical an- 
gular velocity for high fields is completely explained by 
the length-scale analysis done here. However, before ac- 
cepting this result, one needs to check that the small 
imaginary part of the frequency does not significantly al- 
ter it, since the imaginary part of k± is also becoming 
small when the above equation hold s tru e. Consider the 
replacement k — > K-|-i/(r2T^) in Eq. (4.6): the imaginary 
part of K can be ignored as long as 



1 pVi 

> —— 



(5.8) 



i-2 rj ri^ 

Thus, the results presented here are self-consistent for 

'a 



r„ > {%.is)Bl^Tl 



(5.9) 



Finally, now that ds/r) is seen as the relevant thick- 
ness of the VBL when magnetic fields are important, one 
can check that the ohmic dissipation rate within the skin- 
depth region of the crust is smaller than the shear dissipa- 
tion rate within the VBL. Note that the ohmic disspipa- 
tion rate scales as (C/r)|5Jp. For finite conductivity, the 
magnetic field at the boundary is continuous, and thus 
the derivative of the magnetic field in the crust, which 
gives the current, is proportional to {c/47t)SB/(. Thus, 
the ohmic dissipation rate scales as [l/((^r)]|c(5i3/(47r)p, 
and SB ^ [Bk/{Kil.)]Sv. On the other hand, the shear 
dissipation rate in the VBL scales as r]\kSv\'^dB/,y Substi- 
tuting in typical numbers, the ratio of the ohmic to shea r 
dissipation rates is lO^^i^j^^Tio^f^/rio [see Eqs. (| 



(4.12),(4.16)-(4.18)], which justifies the neglect of ohmic 
dissipation in this paper. 
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VI. CONCLUSIONS 

While perhaps not obvious at first, it should now be 
clear why the results found in this paper hold true. It 
has been shown that magnetic fields increase the thick- 
ness of the VBL, which would tend to increase the VBL 
damping time. However, vibrations of the magnetic field 
lines, corresponding to Alfven waves, produce a larger 
shear than one would expect due to the VBL thickness 
alone. The wavelength of these vibrations is the distance 
the Alfven waves travel in one period of the r-mode os- 
cillation. This is a longer distance than viscous diffusion 
can travel in the same amount of time. Thus, the shear 
is smaller than when the magnetic field is zero, but be- 
cause the VBL thickness is increased, the shear acts over 
a larger volume. The net result is that magnetic fields 
decrease the VBL damping time. 

Before concluding this paper, some caveats need to be 
made about the validity of the approximations made in 
this paper. First, the assumption that angular deriva- 
tives are small needs some further comment. The last 
figure show the integrand that gives X for = for 
two cases. This integrand determines the VLB damp- 
ing rate. For the model used in the last section, the 
special angles discussed in Sec. LV, in radians, are lo- 
cated at X_ — 1.23 and Xj^ = 1.91. It can be seen 
that the approximation that angular derivatives are small 
breaks down near X±. However, correcting this should 
not greatly affect the results in this paper, since these 
regions make small contributions the VBL damping rate. 
However, future studies could consider correcting for 
the angular behavior near X±. Second, future studies 
could consider cases where the imaginary part of the fre- 
quency influenc es t he structure of the VBL, as when the 
boun d in Eq. (5_^) is violated. Third, as noted after 
Eq. (4.12), the approximation d/Rc <C 1 breaks down 
for B > 10^^ Gauss. Future studies could consider ex- 
tending the results presented here to the case of flelds 
larger than this bound. More importantly, future stud- 
ies should explore the supcrconducting-superfluid case. 
Mendell |2^ has shown that Alfven waves are replaced 
by cyclotron- vortex waves in a superfluid neutron, type II 
superconducting proton, ordinary-fluid electron plasma. 
Since the square of cyclotron- vortex wave velocity scales 
linearly with B, one would speculate that the results for 
superfluid neutron stars would be similar to those pre- 
sented in Eqs. ( |5.6| ) and (5.7) but with r„ oc B~^/'^ and 
r^c/f^o c>c B^/^^, for large B. However, there are com- 
plications with applying the no-slip boundary condition 
when both the neutrons and protons form inviscid quan- 
tum condensates. In this case, in principle, both the 
neutrons and protons can slip at the crust-core boundary. 
The electrons can slip as well, since the crust is conduct- 
ing. Instead, the boundary conditions are determined by 
pinning effects. In the strong pinning limit neither the 
neutrons or the protons can slip at the boundary. The 
pinning of proton vortices between neutron vortices, the 



possible expulsion of proton vortices during the spindown 
phase of a pulsar, and mutual friction effects due to the 
scattering of electrons off quantized neutron and proton 
vortices, are further effects that need to be considered. 
Thus, future work on the correct boundary conditions is 
needed before the superfluid case can be worked out. 

In conclusion, this paper has shown that magnetic 
fields profoundly affect the VBL in neutron stars. Note 
that this paper treated the crust as perfectly rigid, while 
crustal motion increases the VBL damping rate by a fac- 
tor between 1 and 400 |^9|. Thus, one should treat the 
results presented here as upper bounds on the critical an- 
gular velocity. (However, other effects not yet considered, 
like turbulance, could again increase the critical angular 
velocity.) Still, it is clear that magnetic fields of the mag- 
nitude typically found in neutron stars can substantially 
shorten the VBL damping time and, if sufficiently large, 
can suppress the r-mode instability. 
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